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Abstract 

Much recent progress has been achieved for stabilization of linear and nonlinear systems with input delays that are long and 
dependent on either time or the plant state — provided the dependence is known. In this paper we consider the delay variations 
as unknown and study robustness of nominal constant-delay predictor feedbacks under delay variations that depend on time 
and the state. We show that when the delay perturbation and its rate have sufficiently small magnitude, the local asymptotic 
stability of the closed-loop system, under the nominal predictor-based design, is preserved. For the special case of linear systems, 
and under only time- varying delay perturbations, we prove robustness of global exponential stability of the predictor feedback 
when the delay perturbation and its rate are small in any one of four different metrics. We present two examples, one that is 
concerned with the control of a DC motor through a network and one of a bilateral teleoperation between two robotic systems. 
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1 Introduction 

Networked control systems are present in various en- 
gineering applications such as tele-robotics, remote 
surgery and automotive systems among others. One of 
the major reasons is that they are advantageous over 
traditional control systems in terms of flexibility, relia- 
bility, maintenance cost etc. [14]. However, often their 
performance can be severely degraded in the presence 
of delays induced by the network [3], [13]. When the 
delay is constant and known, the predictor-based con- 
troller compensates the network- induced delay [25]. 
Yet, the networked-induced delay might be subject to 
time-varying and state-dependent uncertainties, which, 
when they are not considered in the control design, 
not only can degrade the performance of the control 
system, but can also destabilize the network [32]. It is 
therefore crucial to quantify the robustness properties 
of the constant-delay predictor feedback in the presence 
of time- and state-dependent delay uncertainties. 

Numerous methodologies exist, dealing with the stabil- 
ity or stabilization of nonlinear systems with input [20] , 
[33], [34], [35], [36], [37], [38], [39], [48], [49], [50], or state 
[16], [18], [21], [22], [23], [42], [45] delays. Predictor-based 
techniques are developed for the compensation of long 
actuator delays in linear [2], [7], [17], [19], [29], [41] or 
nonlinear [4], [5], [24], [31] systems. Among them, [4], 
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[30] , [41] arc dealing with time- varying and [5] with state- 
dependent input delays. Predictor feedback has been 
also successful in designing stabilizing controllers for lin- 
ear systems with uncertainties either on the delays [6], 
[8] or on the plant parameters [40], [47] or on both [10]. 

Although, some of the first predictor-based designs for 
linear, unstable plants, with constant input delays go 
back to the 1980s [2], a Lyapunov construction has been 
unavailable until recently [29]. In addition, a Lyapunov 
functional is provided in [31], which is employed to the 
stability analysis of the proposed control design. Yet, 
the robustness properties of the nominal, constant- 
delay, predictor-based feedback under time- and state- 
dependent delay perturbations remain unexplored. 

We consider forward-complete nonlinear systems that 
are locally, exponentially stabilizable in the absence of 
the delay (by a possibly time- varying control law), for 
which we employ the predictor-based design. The pre- 
dictor controller is designed assuming constant input 
delay and using only an estimation of the unmeasured 
(since the delay is not known) infinite-dimensional ac- 
tuator state. We prove robustness of the constant-delay 
predictor-based feedback, under simultaneous time- 
varying and state-dependent perturbations on the delay. 
Specifically, using the nonlinear infinite-dimensional 
backstepping transformation, we construct a Lyapunov 
functional for the closed-loop system that is comprised 
of the plant, the predictor feedback and the observer 
for the actuator state. With the constructed Lyapunov 
functional, we prove that the closed-loop system remains 
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locally asymptotically stable when the perturbation 
and its rate are small (Section 2). 

We also deal with linear systems under time- varying de- 
lay perturbations. We show robustness of global expo- 
nential stability of the predictor feedback for the cases 
where the delay perturbation and its rate either have 
small magnitude, or their C\ norm is small, or they con- 
verge to zero as the time goes to infinity or, finally, have 
a small moving average for large times (Section 3). Fi- 
nally, we illustrate the robustness properties of the pre- 
dictor feedback with two examples. The first one is an 
example of a DC motor which is controlled through a net- 
work. The network induces a delay which is comprised 
of a known constant part and an unknown time- varying 
perturbation on this nominal value. In addition, the de- 
lay is subject to a state-dependent perturbation that de- 
pends on the armature current. The second example is 
concerned with the bilateral teleoperation between two 
robotic systems through a network. The network induces 
a constant nominal delay which is subject to an unknown 
time-varying perturbation that has a small moving av- 
erage after a long period of time (Section 4) . 

Notation: We use the common definition of class JC 1 /Coo 
and ICC functions from [26]. For an n- vector, the norm 
| • | denotes the usual Euclidean norm. We say that a 
function ( : R + x (0,1) ^ R+ belongs to class K.C if it is 
of class /C with respect to its first argument for each value 
of its second argument and continuous with respect to its 
second argument. It belongs to class AXoo if it is in AX 
and also in /Coo with respect to its first argument. With 
Q we denote the inverse of the function ( with respect to 
its first argument for each value of its second argument. 



2 Robustness to time- and state-dependent de- 
lay perturbations for nonlinear systems 

We consider nonlinear plants of the form 

X(t) = f(x(t),u(t-D-6 (t, X(t)j) ) , (1) 

where / : C 2 (R n x R;R n ), satisfies /(0,0) = 0, D > 0, 
under the nominal, constant-delay predictor feedback 
given by 



u(t) = «(* + A-P(t)) , 



where for all t - D < 6 < t 



P(9)=X(t) + J j[P{s),U{s))ds, 



(2) 



(3) 
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Fig. 1. Control over a network, with delay that varies with 
time (as a result of other users's activities) and may be 
state-dependent. The designer only knows a nominal, con- 
stant delay value D. The delay fluctuation 8(t,X) is un- 
known. 

is the estimate of the actual predictor state P*(9), de- 
fined for aUt-D- 5(t, X(t)) < 6 < t as 



P*{6)=X{t)+ [ 

Jt-D-6(t,X(t)) 

G(a)=5 t (o-(s),P*(s)) 



f(P*(a),U( a )) 
l-G(s) 



+V8(o-(s),P*(s))f(P*(s),U(s))- 



ds (4) 



(5) 



where a, the actual predicted time (which should be 
compared with the estimated predicted time t + D) is 



a(0)=t + 



r 

Jt-D-6(t,X(t)) 1 



ds 



G( s y 



(6) 



An example of such a model together with the predictor- 
based controller is shown in Fig. 1. 

Let us now make clear why P* is the actual predictor 
state of X. Define the actual delayed time as 



(j}{t)=t-D-5{t,X{t)), 

and the actual predicted time as 

a{t)=r\t) 

= t + D + S{a{t),X(a(t))). 



(7) 



(8) 



Then we show that the signal in (4) satisfies P*(t) — 
X(a(t)). Differentiating (8) we get that 



&{t) = l + 5 t (a{t),X(a(t)))&(t) 

+V8(a(t),X(a(t)))X'(a(t))a(t), 



using (1) we have 



and since a(t) - 

F(t) = S t (a(t),X(a(t))) 

+VS(a(t),X(a(t)))f(X(a(t)),U(t))- 



(9) 

(10) 
(11) 
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Define the change of variables t — a (6), where <j>(t) < 
9 < t. With the help of (10) re- write (1) in terms of 9 as 



dP*(9) _ f(P*(6),U(6)) 



d9 



G{9) 



for all <j>(t) <9<t, (12) 



where we substitute X(a(9)) with P*{9). Integrating 
(12) from 4>(t) to 9 we get (4) by using the fact that 
P* (</>(*)) = X(a(<f>(t))) = X(t). Noting that a((j)(t)) = t 
we also get (6). A more detailed discussion about defi- 
nition (4) can be found in [5] . 

The predictor state (3) is the certainty equivalent pre- 
dictor for system (1). This becomes clear by setting 
5 = in (4). Note that P(9) (or P*{9)) should be viewed 
as the output of an operator, parametrized by t, act- 
ing on P(s) and U(s), t - D < s < 9 (or t - D - 
6(t,X(t)) < s < 9) in the same way that the solution 
X(t) to an ODE can be viewed as the output of an oper- 
ator, parametrized by to , acting on X(s) and the input 
U(s), t < s < t. However, P is given implicitly since the 
plant is nonlinear (for the same reason that the solution 
X{t) to a nonlinear ODE is given implicitly). In the case 

of a linear plant X(t) = AX(t) + BU (t-Dj, equa- 
tion (3) for the predictor state can be solved explicitly 
as P{9) = e A ( e - t+£) )x(t) + e^ 6 ^ BU(s)ds, and 

hence, P(t) = e Ab X{t) + e A ^- ^ BU(9)d9. This is 
the standard predictor (used for example in [2]), which 
is obtained using the variations of constants formula for 
the linear ODE satisfied by the plant. An equivalent rep- 
resentation of the signal P(9) is 



p{x, t) = X(t) + D f (p(y, t),u{y, t)) , 



(13) 



for all x G [0, 1], where it, is the estimation of the actua- 
tor state U (9), t — D — 5 (t, X{t)) < 9 < t, which satisfies 



Du t (x,t) =u x (x,t) 
u(l,t) = U(t), 



(14) 
(15) 



that is, 

u(x,t) = U{t + D(x-l)), for all x € [0,1]. (16) 

With this definition, p(x, t) is the output of an operator, 
parametrized by t, that acts on p(y,t) and u(y,t), y e 
[0,x]. With this representation p(l, t) = P{t). 

Note also that from relation (4) we see that for P* to 
be well-defined the denominator in (4) must satisfy the 
following condition for all 9 > t — D (X(t )) 



c>S t (a(9),P*(9)) 

+V5(a(9),P*(6))f(P*(6),U(0)), 



for c G (0, 1], which is a condition on the perturbation 
6, the initial conditions and the solutions of the system. 
As it turns out later on, this condition is satisfied by 
appropriately restricting the perturbation S and the ini- 
tial conditions of the plant. We proceed by imposing the 
following assumptions on the delay-free plant. 

Assumption 1 The plant X = f (A, w) is strongly for- 
ward complete, that is, there exist a smooth positive def- 
inite function R and class K-oo functions a\, ot2 and 
such that 



OR (A) 
dX 



a 1 (\X\)<R(X)<a 2 (\X\) 
f(X,u)<R(X) + a 3 (\u\), 



(18) 
(19) 



for all lei" and for all lo e 



Assumption 1 guarantees that for every initial condition 
and every measurable locally essentially bounded input 
signal, the corresponding solution of the system exists 
for all times. Forward-completeness is a natural require- 
ment for nonlinear plants with input delay. In the ab- 
sence of this assumption, i.e., when the plant exhibits 
a finite-escape time, the control signal might reach the 
plant "too late" . The difference with standard forward- 
completeness from [1] lies in the fact that R in Assump- 
tion 1 is positive definite, in accordance to the fact that 
/(0,0) = 0. 

Assumption 2 There exist positive constants ji, r* , b, 
X* , afunctiona which belongs to class ICoo and a function 
k : C' 3 ([to, oo) xI";R) satisfying for all t >t 



d l+j n(t,C) 



dHd^ 



Ja(KI), o<*<3,;=o 1 

\ f i + am,0<i<3,j = 1...3-i 



such that for the plant X(t) = f (X(t), n(t, X(t))) the 
following holds for all A (to) € D r * , 

\X(t)\ < k|A(t )|e- A * (t -* o) , for all t> to, (21) 
where D r « = {X € R"||A| < r*}. 

Theorem 1 Consider the closed-loop system consisting 
of the plant (1), control law (2), (3) and observer (14)- 
(15). Under Assumptions 1 and 2 there exist positive 
constants c\, c** , class JC^ functions fi, a* , a class AXoo 
function (, and a class ICC function (3 such that if the 
perturbation 5 satisfies 



\Ht,Z)\ + \6 t (t,0\ + \V6(t,0\< Cl +jim.. 



(22) 



for all (t, £) £ [to, oo) x 1™, then for all initial conditions 
which satisfy 



(17) n(t )<c**, 



(23) 
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where 



satisfies 



rt ft 
U(t) = \X(t)\+ I a* (\U(9)\)d9 + U(9) 2 d9 



t-D 



t-D 



U(9) 2 d9, 



(24) 



't-D-max{0,5(t,X(t))} 

it holds that 

H(t) < j8 (C (n(t ), min {r*, c, £>}) , i - t„) , (25) 

/or all t > t and some < c < 1. 

We now introduce the backstcpping transformation. 

Lemma 1 Define the backstepping transformation 

w(x, t) = u(x, t) — k (t + Dx, p(x, t)j , (26) 

together with its inverse, 

u(x, t) = w(x, t) + n(t + Dx, p(x, tfj , (27) 
where p is given for all x £ [0, 1] 6 tTH 



(28) 



System (1) together with the control law (2), (3) can be 
written in the following form 

X(t) - / (X(t), k (t, X(t)) + w(0, t) + u(0, t)) (29) 
Dw t (x, t) = w x (x, t) + ri(aj, t)f(t) (30) 
tf(M)=0, (31) 

where 

f(t)=f (p(0, t),u(0, t) + u(0, t)) - / (p(0, t), u(0, t)),(32) 
and r\ is defined in Appendix A. The observer error 
u{x,t) = u{x,t) — u(x,t), (33) 



p(x,t)=X(t) + D f(p(y,t), 

w(y, t) + K ( t + Dy, p(y, t) \ )dy 



1 An important observation at this point is that p in (13) 
and p in (28) are identical. To see this, observe that, through 
the backstepping transformation, both p and p satisfy the 
same ODEs in the spatial variable x, with the same initial 
condition X. However, p is expressed in terms of the original 
variables (X, u) in the direct backstepping transformation, 
whereas p, is expressed in terms of the transformed variables 
(X, w) and is used in the inverse transformation. 



u t (x, t) = ir{x, t)u x (x, t) — ( 1 — Dn(x, t)^j r(x, t) (34) 
u(l,t) = 0, (35) 

and 

u xt {x,t) = ir(x,t)u xx (x,t) + ir x (x,t)u x (x,t) - (1 

-Dn(x,t)j r x (x,t) + DTT x {x 1 t)r(x,t) (36) 
1 



u x (l,t) 



where 



tt(x, t) = 



7T(M) 



D r(l, t), 



l + x(&{t) - 1) 
a{t) - 1 ' 



(37) 



(38) 



the function a is defined in (8), and the function r is 
defined in Appendix A. Furthermore, 



Dw xt (x,t) =w xx {x,t) +r 2 (x,t)f(t) 

w x (i,t) = -n{l,t)f(t), 

and 



(39) 
(40) 



Dw xxt (x,t) t)+r 3 (x,t)f(t) (41) 

w xx (l,t)=-r 2 (l,t)f(t)+r 4 (t)f(t) + f T (t)r 5 (t)f(t) 

-n(l,t)fp(t)f(p(0,t),u(0,t)+u(0,t)) 

-n{l,t)Dr{0,t)U + Dri{l,t)r{0,t) (1 
df(p(0,t),u(0,t)+u(0,t)) 



-Dw(0,t) 



Oil 



where 



hit) = 



hit) 



-D ri (l,t)7r(0,t)u x (0,t) 
df (p(0,t),u(0,t) + u(0,t)) 
du 



Of (p(0,t),u(0,t) + u(0,t)) 

dp 

df (p(0,t),u(0,t)) 

dp 

df (p(0,t),u(0,t) + u(0,t)) 
du 

df(p(o,t),u(o,t)) 

du 



(42) 



(43) 



(44) 



and r 2 , r 3 , r^ and r§ are defined in Appendix A. 
Proof We re-write system (1) in the form 

X(t) = f(X(t),u(Q,t)) 
u t (x,t) =%(x,t)u x (x,t) 
u(l,t) = U(t), 



(45) 
(46) 
(47) 
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where the actual actuator state u(x,t) (for which an 
observer is given in (14)-(15)) satisfies (46)-(47) and 

u(x, t)^U{<j){t + x {a(t) - t))) , for all x G [0, 1],(48) 

where cf) is defined in (7), or, incorporating 5, as 

u{x, t) = u(t + x(b + 5 (a(t), X(a(t)))^j - D 

-8(t + x(p + 6 (a(t),X(a(t)))) ) ) . (49) 

The rest of the lemma is proved with lengthy but 
straightforward calculations and using Lemma 11 from 
Appendix B, the backstepping transformation (26), its 
inverse (27), and (13), (14)-(15), (45)-(47), (28), (33). 

For 7r to be a meaningful propagation speed it should be 
positive and uniformly bounded from below and above. 
Using (38) , one can conclude that 5, the initial conditions 
and the solutions of the system should satisfy, in addition 
to (17) which guarantees that < &(t), also 

< D + 6(a(6), P* (9)), for all 6 > t - D (X(i )),(50) 

which guarantees < a(t) — t. The two conditions (17), 
(50) incorporate the functions a and P* , that is, they 
are not expressed in terms of the perturbation <5 and the 
functional II. We derive next a sufficient condition for 
(17), (50) to be satisfied, in terms of II. 

Lemma 2 There exist positive constants c\, c* , such 
that if the perturbation S satisfies (22), then for all solu- 
tions of the system satisfying, 

U{t)<c*, (51) 
they also satisfy 

( Cl + A (\P*(0)\)) (1 + |/ (P*(9), U(6))\) < R, (52) 
for all <j)(t) < 9 <t, where 

i? = min{r*,c,£>}, (53) 

for some < c < 1, and hence, conditions (17) for 
< c < 1, and (50) are satisfied. 

Proof See Appendix C. 

Lemma 3 There exist a continuously differentiable, 
positive definite function S* , a class /Coo function a* 
and positive constants X, c\ , c* such that if the perturba- 
tion 6 satisfies (22) then for all solutions of the system 



satisfying ( 51 ), the Lyapunov function 

V(t)=S*(X(t))+g n f e 9lX \u(x,t)\dx + g 6 f e^ 
Jo Jo 

xu x (x,t) 2 dx + g 9 D / e 9wX \w x (x,t)\ dx 
Jo 

+g 12 D [ e 93X a*(\w{x,t)\)dx + g 8 D f e 9 * x 
Jo Jo 

xw xx (x,t) 2 dx + g 7 D e 9iX w x (x,t) 2 dx, (54) 
Jo 

where gi > 0, i = 1 . . . 12, satisfies 

V(t) < y(f )e~ A(t ~ to) , for all t > t . (55) 
Proof See Appendix D. 

The next two lemmas relate the Lyapunov function V 
with the norm of the system in the original variables, 
represented with PDEs, and the norm in PDE represen- 
tation with the norm in standard delay form. 

Lemma 4 There exist a positive constant c* , a class 
/CCoo function a 2 A and a class /Coo function a 2 5 such that 
for all solutions of the system satisfying (51) the follow- 
ing holds 



«24 (T(t),R) < V(t) <a 25 (T(t)), 
where 



(56) 



T(t) = \X(t)\+ f a* (\u(x,t)\)dx+ f u x (x,t) 2 dx 
Jo Jo 



+ / u x (x,t) 2 dx + / u xx (x,t) 2 dx. 
Jo Jo 

Proof See Appendix E. 



(57) 



Lemma 5 There exists positive constants c* , C\ and 
class /CCoo functions £i and ( 2 such that if the perturba- 
tion 5 satisfies (22), then for all solutions of the system 
satisfying (51) the following holds 

( 1 (n(t),R)<T(t)<( 2 (Il(t),R). (58) 
Proof Using Lemma 2 and (8), (10) we get that 



1 1 1 
< — < 



(59) 



D + R <r{t) -t D-R 

< v{0) < t—^> for a11 <H*) < 9 < L (6°) 



1 + R 



1-R 1 



With relations (16), (48), (24), (57) and applying the 
appropriate change of variables in the integrals the proof 
is immediate using (59), (60). 
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Proof of Theorem 1: Using Lemma 5 we conclude that 
(51) is satisfied if T(t) < d(c*,R), and hence, with 
Lemma 4, (51) is satisfied if 

V(t)<a 24 (Ci(c*,R),R), (61) 

is satisfied. Assume for the moment that (61) is satisfied. 
With Lemmas 3, 4 and 5, relation (61) is satisfied if c** 
in (23) is such that 

c** < ( 2 Ks 1 (a 24 (Ci (c* ,R),R)),R). (62) 

Using (55), with some routine class IC majorizations that 
involve Lemmas 4 and 5, we get estimate (25). 

3 Robustness to time-varying delay perturba- 
tion for linear systems 



is a consequence of the fact that when S does not de- 
pend on the state, the conditions (17), (50) are satisfied 
without restricting the supremum norm of the real ac- 
tuator state U(0), for all 4>(t) < < t. Secondly, in the 
linear case the control law, as well as the direct and in- 
verse backstepping transformations are given explicitly 
and are globally well-defined. 

When (65) does not hold, one can still derive exponential 
stability of the closed-loop system, by imposing other 
conditions on the perturbation S, such us the ones from 
[26] Chapter 9.3. However, in this case one has to guar- 
antee in addition that the propagation speed tt is still 
uniformly bounded from above and below and strictly 
positive. Therefore, we make the following assumptions 
which S has to a priori satisfy. 

Assumption 3 The perturbation S satisfies 



We consider the following special case of system (1) S'(t) < 1, for all t>t , (68) 



X(t) = AX(t) +BJj(t-D- 5{tf) , (63) and is such that 

where for the rest of the section S is a function only 7r* = — — S'(9)) > ^ 

of the time t. For this linear case, the predictor-based su Pe>o-(t )l v I) 

controller is given explicitly as 

Assumption 4 The perturbation 5 satisfies 
U(t) = Ke^X(t) + £_^*-»BU{O)d0. (64) f) + S( t )>0, forallt>t Q , (70) 



Theorem 2 Consider the closed-loop system consisting 
of the plant (63), observer (14)-(15), and control law 
(64)- There exists a positive 5\, such that if the pertur- 
bation S satisfies 



\5(t)\ + \5'(t)\<5 u for all t> t , 



(65) 



then, the closed-loop system is exponentially stable, in 
the sense that there exist positive constants R and A such 
that the following holds: 



n L (t) < RTl h (t )e 



-A(t-to) 



n L (t) = \x 



(t)\ 2 + I 

Jt-D- 



max{0,<5(*)} 



U(9) 2 d6 



U(9) 2 d9. 



(66) 



(67) 



It-D 



id is such that 



™Ve>*(t )(D + 5(9)) 



> 0. 



(71) 



Theorem 3 Assume that 5 satisfy Assumptions 3, 4- 
There exist positive 82 and S3 such that if the perturbation 
S satisfies either 



r {\5'(o)\ 

Jo 



+ \6{9)\)d9<5 2 , 



or 



\6(t)\ + \5'(i)\->0, whent^f 00, 



(72) 



(73) 



or 



The proof of Theorem 2 is based on the application of 
Lemmas 1, 3, 4 and 5 for the special case of plant (63) 
and the special perturbation <5 as in (63). However, we 
give each of these lemmas specialized to the present case 
for two reasons. Firstly, in the special case of linear sys- 
tems with only time- varying perturbation, we study sta- 
bility of the closed-loop system in the H\ norm. This 



1 r t+A 

-^J (\5'(9)\ + \S(9)\)d9<S 3 forallt>T, (74) 

for some positive A and nonnegative T, then, the closed- 
loop system consisting of the plant (63), observer (14)- 
(15), and control law (64) is exponentially stable, in the 
sense that there exist positive constants R and A such 
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that the following holds: 

n L (t)<i?n L (t )e- A(t -* o) 



n L (t) = |x(t)r 



t-£>-max{0,i5(t)} 



U{0) 2 d6 



U(9) 2 d0. 



(75) 



(76) 



From Theorem 3 we observe that in the case of linear sys- 
tems and when the perturbation depends only on time 
and not on the state, we prove robustness of global ex- 
ponential stability of the predictor feedback under three 
alternative conditions on the delay perturbation rather 
than just restricting the magnitude of the delay and its 
rate. This is not possible in the case where the perturba- 
tion depends on the state because for the conditions (17) 
and (50) to be satisfied one has to necessarily restrict 
the magnitude of the perturbation and its rate such that 
they are both sufficiently small. 

We introduce now the backstepping transformation of 
the estimated actuator state. 

Lemma 6 Consider the backstepping transformation 

w(x,t) =u(x,t) - Ke Abx X(t) 

—DK [ e Ab{x -^Bu{y,t)dy, (77) 
Jo 

together with its inverse given by 

u(x, t) = w{x, t) + Ke {A+BK)bx X{t) 

+DK [ X e^ A+BK ^ x -^Bw(y,t)dy. (78) 
Jo 

System (63) together with the control law (64) can be 
represented as 

X(t) = (A + BK) X(t) + Bw(0, t) + Bu(0, t) (79) 

Dw t (x, t) = w x (x, t) - DKe Abx Bii(0, t) (80) 
u>(l,t)=0, (81) 



where the observer error 
u(x,t) — u(x,t) — u(x,t), 
satisfies 



(82) 



Ut(x,t) =Tr(x,t)u x (x,t) — ^1 — bi:(x,t)j r(x,t) (83) 
u(l,i)=0, (84) 



with 



w(x, t) 



l + x(a(t)-l) 
a(t) - t 



4>(t) =t-D- S(t) 

o-{t)=r\t) 

= t + D + 5(a(t)) 
1 



(86) 
(87) 

(88) 



r(x, t) = -jrw x {x, t) + Ke {A+BK)Dx (A + BK)X(t) 

+KBw(x, t) + DK(A + BK) 

x f e^ A+BK ^ x -y^Bw{y,t)dy. 
Jo 

Furthermore, 

Dw xt {x, t) = w xx (x, t) - b 2 Ke ABx ABu{0, t) (89) 
w x {l,t) = DKe Ab Bu{Q,t). (90) 

Proof System (63) can be re-written in the form 

X(t)=AX(t) + Bu(0,t) (91) 
u t (x,t) =n(x,t)u x (x,t) (92) 
u(l,t) = U(t), (93) 

where tt is defined in (85). The predictor feedback (64) 
can be written as 

U(t) = Ke Ab X{t) + D [ e Ab ^-y^Bu{y, t)dy, (94) 
Jo 

where the estimation of the unmeasured actuator state 
U(0), for alU - D - S(t) < 6 < t, u(x,t) is defined in 
(14), (15) and satisfies (16). With representation (91)- 
(93) for system (63), the actuator state u(x,t) is 

u(x, t) = U((j)(t + x (a(t) - t))) , for all x e [0, 1],(95) 
or, incorporating S, u(x,t) = U 

-t) - 5 (t + x (t) + £(cr(i))))) for alia; G [0,1]. Since 

the perturbation 5 satisfies (65), it follows from defini- 
tion (85) and relations (86), (87) that o-(t) — t > and 
that 1 — S'(t) > 0, for all t>t n . Define the quantities 



(96) 
(97) 



sup e >„ (to) (l-5'(0)) 

7Tn = ~ . 

sup e > a(to) (D + 5(0)) 



(85) 



From (65) it follows that sup 0>(T ( O )(l — S'(8)) < oo and 
that sup e>a ^(b + 6(9)) < oo, and hence, tt" > 0, 
TT^* > 0. Since a(t) - t = D + S(a(t)) and &(t) = 
i-s\a(t)) ' usm S (65) we conclude that 7r is positive and 
uniformly bounded from above and below. Hence, 7r is a 
meaningful propagation speed. The rest of the proof is 
based on algebraic manipulations and it is omitted. 
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Lemma 7 There exist positive constants n and r 2 such 
that the derivative of the Lyapunov function 

V L {t) = X{t) T PX(t) + h [ e bx u(x, tfdx 

Jo 

+Db 2 / (1 + x)w(x,t) 2 dx 
Jo 

+Db 2 (1 + x)w x (x,tfdx, (98) 



along the solutions of (79)-(81), (83)-(84), (89)-(90) 
satisfies 

Vh(t) < -riV h (t) + r 2l (t)V L (t), (99) 
7 (t)=max{|5(a(t))|, 

6(a(t))(l-S'(a(t)))-D6'(a(t))\}. (100) 
Proof See Appendix F. 

Lemma 8 There exists a positive S\ such that if the per- 
turbation S satisfies (65), then there exists a positive A 
such that V in (98) satisfies 



V h (t) < -\V h (t). 
Proof See Appendix G 



(101) 



Lemma 9 There exist positive constants M 5 ^ andM§^ 
such that 

M 5iL V L (t) < T h (t) < M 6;L Vl(*) (102) 

r L (i) = \X(t)\ 2 + / u(x,t) 2 dx + u(x,t) 2 dx 
Jo Jo 

+ / u x (x,t) 2 dx. (103) 
Jo 

Proof See Appendix H. 

Lemma 10 There exist positive constants M-j^ and 
Mg^L such that 

M 7 , L r L (i) < n L (t) < M 8;L r L (i). (104) 

Proof From (16) we get u x (x,t) = DU'(t+D(x - 1)). 
Applying a change of variables in (103), with (95) we get 

rt 



r L (i) = |*(*)r + 



1 



i r 

W) - t Jt-D 



i 



t-D-6(t) 
t 



-U(6) 2 d6 



7: / U(6) 2 d6 + 

D Jt-D 



[ U(9) 2 d6. (105) 

Jt-D 



Hence, the lemma is proved with M-j^ 
M L 



max{l,M L + A}' 

SU Pt-b-S(t)<0<t l-sV(fl)) yr _ 1 

inf t > t0 (Zj+«(<r(t))) ' 8 < L ~ min{l, ff 5* Wl ".^}' 



Proof of Theorem 2: Using Lemma 8 and the comparison 
principle ([26]), we get V h (t) < y L (t ) e - A ( t ^ t °). With 

Lemmas 9 and 10 we get (66) with R = 



M 6iL M 8 , L 

M 5 , L M7,L ' 



□ 



Proof of Theorem 3: We consider first the case where S 
satisfies (72). Under Assumptions 3 and 4, Lemmas 6, 
7, 9 and 10 apply to this case as well (with 7Tq* and 7r** 
replaced by ttq and tt'{ respectively). The only difference 
with the proof of Theorem 2 is in the proof of Lemma 8. 
Towards that end, we solve (99) to get 



V L (t) < e J *o VL(*o) 



(106) 



Consider that j(t) 



7r(0,t) 



D 



Applying the change 



of variables r = (j>(9) in the integral and using the facts 
that a (t ) > t and <j>'(t) = 1 - S'(t) we get 

/ l(r)dr<- \5(6)\d6<^. (107) 

Jt ^1 J to n l 



Analogously, for the case 7(f) 



[ 



-f(r)dT < 



D 



we get 
(108) 



Hence, Lemma 8 is proved with r 2 S 2 < min{7T*ri , 
and the fact that f ' ^(r)dT < f t °° j(r)dr. 



n+K*~D) [ """" J t 'V/"' — jt 

Note that in the present case, Lemma 8 can be proved 
directly from relation (99) using Lemma B.6 in [28]. 

For the case where 5 satisfies (73), Lemma 8 is proved 
by combining Lemma 7, Lemma B.8 in [28] and the fact 
that 7(i) satisfies 



1 



7Tl 



j(t)<(l + -+D)(\S(a(t))\ + \S'(a(t))\) 



(109) 



Finally, if S satisfies (74), Lemma 8 is proved using 
Lemma 9.5 in [26]. □ 

4 Examples 

4-1 Control of a DC motor over a network 

We consider the following model of a field-controlled DC 
motor ([44]) with negligible shaft damping 



duj(t) 

dt 
di a (t) 

dt 
di f (t) 

dt 



= 6k(t)i a (t) (110) 

=-bi a ,(t) + k-ci i (t)u(i) (111) 
=-aii(t) + U {t-D- p(t,i f (t)Mt)^{t)))^) 
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where if, i a arc field and armature currents respec- 
tively, u> is angular velocity and a, 6, c, are positive 
constants. The equilibria of the unforced system arc 
(w,i a ,if) = (wo,f,0) for some constant cjq- The sys- 
tem is feedback linearizable for (w,i a ,if) e D, D = 
{ (w, i a , if ) € M 3 |cj > and i a > ^} . A delay- free design, 
based on full-state linearization, is (Chapter 13.3 in [26]) 

U(t) = - (-JfiZi(t) - K 2 Z 2 {t) - K 3 Z 3 {t) - a) (113) 

7 

Zi(t) = #z a (t) 2 + cw(t) 2 -0^- clj 2 (114) 

Z 2 {t)=20i a (t)(k-bi a {t)) (115) 
Z 3 (t) = 29 {k- 2W a (i)) (-W B (t) + fc - ci f (t)w(t)) (116) 
7 = -2c0(fc- 26i a (i))w(t) (117) 
a = 2ca6> (fc - 26i a (f)) i f (t)w(t) - 260 (3fc - 46i a (t) 
-2ci f (t)w(t)) {-bi a (t) +k- ci { (t)uj(t)) 
-2c9 (k - 2bi a (t)) if(t) 2 uj(t). (118) 

Shifting the equilibrium (w 0l f > 0) of the system to the 
origin and setting Xl = w - lj , X 2 = i a - |, X 3 = i f , 
6(t,X 1 {t),X 2 (t),X 3 (t)) = p(t,i { (t),i a (t),oj(t)) we get 

at 

Xl(t)=flX 2 (t)X 3 (t) + yX 3 (t) (119) 

X 2 (t) = -bX 2 (t) - cX^XAt) - cu X 3 (t) (120) 
X 3 (t) = -aX 3 (t) 

+u(t-D-S(t,X 1 (t),X 2 (t),X 3 (t))) . (121) 

The motor is controlled through a network that induces 
a constant delay D (e.g. [11]). The known, constant de- 
lay, is subject to a time- varying perturbation due to the 
effect of transmission of control signals to other motors 
through the network. We further assume that the pertur- 
bation S increases when the armature current increases. 
Define the estimated predictors of X\ , X 2 and X 3 as 

P 1 (t)=X 1 (t) + eJ ^ (p 2 (s)P 3 (s) + ^P 3 (s)^ da(122) 

h(t)=X 2 (t) + J (-bh(s) - cA (*)£»(*) 

-cuj P 3 (s)^ds (123) 
P 3 {t)=X 3 (t) + (-aP 3 {s) + U(s)^ds, (124) 

respectively. Setting in (113)— (118) u> = X\ + ujq, i a = 
X 2 + 1 , if = X 3 and replacing X\ , X 2 , X 3 by the predic- 
tors (122)-(124) we get the nominal predictor feedback. 

We choose the set-point for the angular velocity of the 
motor as ojo — 2, the nominal delay D = 1 and the pa- 



0.3 




Fig. 2. The field (top) and armature (bottom) currents for 
a network controlled DC motor modeled by (110)— (112) 
with the nominal predictor feedback, under an input de- 
lay perturbation 5(i,i a (t)) — 0.5i a (i) 2 + 0.2sin(£) 2 (solid 
line) and <5(t,i a (i)) = (dotted line). The initial condi- 
tions are i f (0) = 0.1, i a (0) = 0.8, w(0) = 1 and U(6) = 0, 
-l-5(0,i a (0)) <6»<0. 

rameters of the plant asa = b = c = k = 9 = l. The 
delay perturbation is S(t,X 2 (t)) = 0.5 (X 2 (t) + \ ) 2 + 
0.2sin(i) 2 . The initial conditions for the plant and the 
actuator state are chosen as Xi(0) = —1, X 2 (0) = —0.2, 

X 3 (0) = 0.1 and U{9) = 0, -1 - 0.5 (X 2 (0) + l) 2 < 9 < 
respectively. The parameters of the controller are cho- 
sen as K\ = — 1, K 2 = K 3 = — 3, such as the lineariz- 
able, delay- free, system (i.e., the delay- free plant in the Z 
coordinates) has three eigenvalues at —1, and the initial 
estimate of the actuator state as U(9) — 0, — 1 < 9 < 0. 

In Fig. 2 we show the field and armature currents, and 
in Fig. 3 the input voltage and the angular velocity of 
the motor. The nominal predictor feedback achieves lo- 
cal stabilization of the closed-loop system at the desired 
equilibrium, despite the presence of the perturbation. 

4-2 Bilateral teleoperation 

In bilateral teleoperation [15], the operator (e.g. a hu- 
man) controls a robotic system, called the master, at the 
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Fig. 3. The angular velocity (top) and the field voltage 
(bottom) for a network controlled DC motor modeled by 
(110)— (112) with the nominal predictor feedback, under an 
input delay perturbation 5(t,i a (t)) = 0.5i a (t) 2 + 0.2sin(t) 2 
(solid line) and S(t, i a (i)) = (dotted line). The initial con- 
ditions are i f (0) = 0.1, i a (0) = 0.8, w(0) = 1 and U(0) = 0, 
-l-<5(0,i a (0)) < 6» < 0. 

one end of the communication network. The actions of 
the master are transmitted (through the network) to an- 
other robotic system, called the slave, at the other end 
of the network. The goal of the control algorithm is the 
slave manipulator to behave (in a certain sense) as the 
master manipulator. A model of two robotic systems, 
each one having n degrees of freedom, representing the 
master and the slave manipulators is written as ([15]) 

x m (t)+x m (t)=r m (t-D-6(t)) (125) 
x s (t)+x s (t)=T s (t-D-26(t)Y (126) 

where x m , x s € R" are the degrees of freedom of the 
robotic systems and the torques r m , r s €E R™ are to be 
designed such as coordination between the master and 
the slave is achieved asymptotically, i.e., x m — x s — »■ 
as t — > oo. The constant delay D represents the known, 
network-induced delay which is subject to time- varying 
perturbations that are often present due to congestion, 
distance etc. [12]. For simplicity we assume scalar x m , 



x s , T m , r s and we re-write (125), (126) as 

X 1 (t)=X 2 (t) (127) 

X 2 (t) = U 1 (t-D-5(tj) (128) 

X 3 (t)=X 4 (t) (129) 

X A {t) = U 2 (t-D- 25(t)) , (130) 

where X\ = x m , X 2 = x m , X 3 = x a , X A = x s , U\ = 



r m and U 2 — t b . A simple controller is ([15]) r m (i) = 
-K p (x m (t) -x s (t)) - B m x m {t) - K p (x m (t) - r) and 
r s (t) = K p (x m (t) - x s (t)) - B s x s (t) - K p (x s (t) - r), 
where r is the set-point for the positions of the manipu- 
lators. The predictor-based version of this controller is 

E/i(i) = ~K P (A(t) - P 3 (t)) - B m P 2 (t) 

-^p(AW-r) (131) 
U 2 {t) = K p (A(t) - P 3 (t)) - B s P 4 (t) 

-K p (P 3 (t)-r) (132) 

Pi(t)=Xi(t)+ [ P t+1 (6)dd, 1 = 1,3 (133) 

Jt-D 

Pj (t) = Xj (t) + f Ui (6)d6, j = 2, 4. (134) 

Jt-D 2 

We choose the desired set-point for x m and Xg as v — ^, 
the parameters of the controller as K p = B m = B s = 2, 
the known delay as D = 1, the initial condition of the 
plant as x m (0) = i m = x s = 0, x s (0) = 1, the initial 
actuator state as 17(0) = 0, -D - (5(0) < < and 
the initial estimation of the actuator state as ?7(0) = 0, 
— D < 9 < 0. We illustrate the robustness properties of 
the predictor feedback under a time-varying delay per- 
turbation that is neither in C\ nor converges to zero as 
the time goes to infinity nor is small in magnitude. How- 
ever, after some long period of time its mean is small. 
This disturbance is show in Fig. 4. 

In Fig. 5 we show the difference between the positions of 
the master and the slave for the cases where either there 
is or there is not a perturbation <5. In both cases, under 
the nominal predictor feedback the position of the slave 
tracks the position of the master. In Fig. 6 we show the 
torques applied to the two robotic systems under the 
perturbation 5. The control efforts are oscillatory as a 
result of the effect of the oscillatory perturbation. 

5 Conclusions 

Looking into the details of the proofs, we note that in the 
case of nonlinear systems with state-dependent pertur- 
bations, there is a trade off between the achievable re- 
gion of attraction and the size of the perturbation and its 
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The delay perturbation S satisfying 
-S(t) + 0.1sin(t) 2 ,5(0) = 1 induced by the network 
The model of the two robotic 



Fig. 

5(t) 

in bilateral teleoperation. 
systems is (125)-(126). 




Fig. 5. The error between the position of the master and 
the slave manipulators for two robotic systems modeled by 
(125)-(126). The two manipulators are coordinated through 
a network with the predictor feedback (131)-(134), under an 
input delay perturbation S(t) = (dotted) and 8 satisfying 
5{t) = -6{t) + 0.1sin(t) 2 , 5(0) = 1 (solid), induced by the 
network. The initial conditions are a; m (0) = 0, x a (0) = 1, 
±m(0) = i B (0) = 0andr m (6>) = r s (6>) = 0, -1-8(0) < 9 < 0. 

rate, at the origin. For linear systems under time- varying 
perturbations, global exponential stability holds, but the 
size of the perturbation and its rate should be appropri- 
ately restricted. With the available Lyapunov functional, 
our next step is to study the inverse optimal redesign 
problem of predictor feedback for nonlinear systems. 

One might raise the question of robustness to stochastic 
delay perturbations, since stochastic perturbations have 
some resemblances with the time-varying case. Yet, in 
our analysis we restrict not only the magnitude of the 
perturbation 5 but also the magnitude of its derivative 
(which also guarantees the invertibility of <j> = t — D — 5 ) , 




Fig. 6. The input torques of the master and the slave ma- 
nipulators for two robotic systems modeled by (125)-(126) 
coordinated through a network with the predictor feedback 
(131)-(134), under an input delay perturbation 5 satisfy- 
ing S(t) = -S(t) + O.lsin(i) 2 , 5(0) = 1, induced by the 
network. The initial conditions (0) = 0, x B (0) = 1, 

im(0) = ± s (0) = and r m (<9) = t s {6) = 0, -1-5(0) < 6 < 0. 

which can be unbounded in the case where S is white 
noise, or even when S is a low pass version of white noise. 

Appendices 

A The perturbation signals of Lemma 1 

With u defined in (27) in terms of p and w, the pertur- 
bation signals n, r, r 2 , r 3 , r 4 and r§ are 



n{x,t) = -D 



dn (t + Dx, p(x, t)j 

d~ P 

3f(p(i/,t).a(y,t)) . 



xe Jo " (A.l) 

2 8k ( t + Dx, p(x, t) J 
r(x, t) = -j-w x (x, t) H ^ — '- 



+ 



8k (t + Dx, p(x, t)j 
8~ P 1 



{p{x,t),u{x,t)) (A.2) 



r 2 (x,t)=-D D- 



8 2 k 



(t + Dx,p(x,t)j 



dpdt 



+f T (p{x,t),u{x,t)) 



8 2 n(t + Dx,p(x,t) 



dp 



n2 



Xe Jo *t ly _ D 2 

On (t + Dx, p(x, i)) g f (/5(a , 5 t j } ^ t jj 



r 3 (x,t) =r 2x (x,t) 



dp 



dp 



(A.3) 
(A.4) 
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r 4 (t) = D 



d 2 K(t + D,p(l,t)) - fl 
9 2 ^ i-e Jo 



3/(p(x,t),a(x,t)) 



dx 



+D 2 f(p,u) 



dtdp 



d 2 p 



xe Jo a? -| i — '- 



xe Jo 55 D 



dp 



(I 



'- 1 d 2 f{p{x,t),u{x,t)) , , 

<| / ^ f{p{x,t),u)dx 

1 d 2 f(p(x,t),u(x,t)) 



+D 



r(x,t)dx] (A.5) 



r 5 (t) =e 



£> f 1 

— e J o 



<9/S<9m y 

_ , t),H(y,t)) j T 8 2 K (t + D,p(l,t)) 



dp 



xD e Jo ^ 



d 2 p 



(A.6) 



where the notation 9 /(p( ^MjK"^.*)) f(p(x 1 1), u(x, t)) 
corresponds to a matrix Q = {gi,j} 1<:i j <n with ele- 
ments q itj = 9 h{p{ g^ {x,t)) f{p{x, t),u{x, t)), where 
/ = {h, ■ • ■ , fn) T and p = (pi, . . . , p„) T . 



B Technical Lemmas 
Lemma 11 T/ie predictor p in (13) satisfies 
Dp t ( a; ,t)=p x ( a; ,i) + J De Z3 Jo a, * /(t)) (B .i) 
w/iere / is defined in (32). 

Proof Differentiating (13) with respect to t, x and using 
(14)-(15) with the fact that p{0, t) = X(t) we get 

*{x,t) = Df{p{0,t),u{0,t)+u{0,t)) 



+D 



df {p{y,t),u{y,t)) 



du 

-Df{p{x,t),u{x,t)) 
y{x,t) = Dp t {x,t)-p x {x,t). 



u y {y,t)dy 



(B.2) 
(B.3) 



Since 



f{p{x,t),u{x,t)) ■ 



a 



df (p(y,t),u(y,t)) 
dp 



p y {y,t) 



df{p( y ,t)My,t)) , 
h — 



9u 

+/(p(0,t),fi(0,t)), 



(B.4) 



we get that 



*(x,t) = D [' 
Jo 



df {P{y,t),u{y,t)) 
dp 



+ D/(t).(B.5) 



Solving (B.5) for VP the lemma is proved. 

Lemma 12 There exists a class K-oo function ae such 
that for all x € [0, 1] 

\p{x,t)\ <a 6 (\X{t)\+J a*(\u(x,t)\)d x y (B.6) 

Proof Differentiating (13) with respect to x and com- 
paring the resulting ODE with the ODE in t for X, the 
proof is complete with Assumption 1 and the comparison 
principle after appropriately majorizing e D Da^ < a*. 
The detailed proof can be found in [31] (Lemma 7). 

Lemma 13 There exists class /Coo function an . . . a i3 
such that for all x e [0,1] 



\w{x,t)\<a n {n{t)) 
\w x (x,t)\ < \u x {x,t)\ +a 12 



(B.7) 
(B.8) 



/ w xx (x,t) 2 dx <6 u xx (x,t) 2 dx + Qi 3 ,(B.9) 
Jo Jo 



where 



fl{t) = \X{t)\+ [ a*{\u{x,t)\)dx+ [ u x (x,t) 2 dx.(B.10) 
Jo Jo 

Proof The proof of the lemma is based on algebraic ma- 
nipulations and routine class JC majorizations using the 
direct (26) backstepping transformation together with 
relations (13) for the predictor state and Lemma 12. For 
the reader's benefit we prove (B.7) and (B.8). The rest 
can be proved similarly. From (26) and (20) we get that 
\w{x, t)\ < \u{x, t)\+ a {\p{x, t)\). Using the fact that 

sup \u{x,t)\ <\u{l,t)\+ [ \u x (x,t)\dx, (B.ll) 
xe[o,i] Jo 

with relation (2) and Lemma 12 we get (B.7). For proving 
(B.8) we proceed as follows. Differentiating (26) with 
respect to x we get that 

8k It + Dx,p(x, t) \ 
w x {x, t)=ii x {x, t) + D ^ '- + D 



dt 



8k (t + Dx,p{x, t)j 



dp 



■f{p{x,t),u{x,t)).{B.12) 



Combining (20), (B.32) with Lemma 12 and (B.ll) we 
arrive at (B.8) with appropriate class JC majorizations. 
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Lemma 14 There exists positive constants M* , c* such 
that for all solutions of the system satisfying (51) the 
following holds for all x G [0, 1] 



(x,t)\ < M* ( \X 



:(t)\+ fa*( 
Jo 



\w(x,t)\)dx . (B.13) 



Proof Under Assumtpion 2 and choosing c* < R, from 
Theorem 4.14 from [26] there exist a continuously dif- 
ferentiable function £ : [t ,oo) x Dr — > K n , where 
D R = {X e W l \\X\ < R}, and positive constants Mi, 
M 2 , M 3 and M 4 such that for all X e D R 



Mi\X\ 2 <S(t,X) < M 2 \X\ 2 
dS(t,X{t)) 8S(t,X(t)) 



dt 



dS (t,X) 



< 



dX 
-M 3 \X(t)\ 2 



f(X(t),K(t,X(t))) 



dX 



<M 4 \X\. 



(B.14) 

(B.15) 
(B.16) 

for all X e Dft and every 
co E K such that \uj\ < M for some positive constant M, 
there exists an increasing function in both arguments 
L e C M+) such that along the solutions of X(t) = 
f {X(t),K (t, X(t) + uj(t))) it holds that 



Since feC 2 



S<-M 3 \X(t)\ 2 + y gx > (f(X(t), K (t,X(t))+u(t)) 

-f(X(t),n(t,X(t)))) 
<-M 3 \X(t)\ 2 + Mi\X(t)\L(R,M)\Lj{t)\, (B.17) 

where we used Lemma 3.1 in [26]. With S* = VS we get 
^( t ,X(t))<--^|^( t )| + ^S^| W ( t )|.(B.18) 



2VMT 



2VMT 



Differentiating (28) with respect to x we get that 

p x (x,t)=Df (p(x,t),n(t + t)x,p{x,t)^ 
+ w(x,t)) . 



(B.19) 



Using the fact that for all x € [0,1], \w(x, t)\ < 
Jp 1 \w x (x,t)\dx (which follows from (31)), relation (B.8) 
together with (51) and Lemma 5 give that for all 
x e [0,1], \w(x,t)\ < M,withM = ( 2 (c*) + a 12 ((2(c*)). 
Using a change of variables in (B.19) as x' = t + Dx and 
comparing the resulting ODE in x' for p with the ODE 
in t for X(t) = f (X(t),K(t,X(t)+oj(t)), with (B.14) 
and after appropriately majorizing s < a*(s), the proof 

is complete with M*(R, M) = ^ + D Mi ^ M) , and 
hence, with c* = min{_R, c*}, where c* satisfies 



M*(R,M(cl))(cl+ ail (( 2 (cl)))<R. 



(B.20) 



Lemma 15 There exists class JCCoo functions an . . . aie 
and a positive constant c* such that for all solutions of 
the systems satisfying ( 51 ), the following holds 

\u{x,t)\<a u (Y{t),R) (B.21) 
\u x (x,t)\ < \w x (x,t)\+a 15 (Y(t),R) (B.22) 

t) 2 <6 t) 2 dx+a 16 (F(t),i?),(B.23) 

Jo Jo 

for all x e [0, 1], where 

Y(t)=\X(t)\+ f a*(\w{x,t)\)dx+ [ w x {x,t) 2 dx.{B.24) 
Jo Jo 

Proof Choose c* as in Lemma 14. Then, the proof of the 
lemma is based on algebraic manipulations and routine 
class JC majorizations using the inverse transformation 
(27), relation (28) for the predictor state and Lemma 14. 

Lemma 16 There exist class AXco functions an . . . a 23 
and positive constants c* , pi , p 2 , p 3 , /14 , /15 such that for 
all solutions of the system satisfying (51) the following 
holds for all x € [0, 1] 



/' 

Jo 

f 

Jo 



\r(x,t)\ < jr\w x (x,t)\+a 17 (Y(t),R) (B.25) 

(B.26) 
(B.27) 

(B.28) 

(B.29) 
(B.30) 



|ri(jc,t)|</ii+ai8 (Y(t),R) 
\r 2 (x,t)\<p 2 +a lg (Y(t),R) 

r 3 (x, tfdx <p 3 + a 20 (Y(t), R) 

\r 4 (t)\<p 4 + a 21 (Y(t),R) 
\r 5 (t)\ < p 5 + a 22 (Y(t),R) 

r x {x,t) 2 dx<a 23 (Y(t),R) 



+ 7v7 / w xx (x,t) 2 dx, 

LJ 2 - Jo 



(B.31) 



where Y(t) is defined in (B.24)- 

Proof Let c* be as in Lemma 14. The proof is based 
on (A.1)-(A.6) combined with (20), the fact that / is 
twice differentiable and with similar calculations as in 
the proof of Lemma 13. Yet, we provide the proofs of 
(B.25), (B.31) as some of the steps are useful later on. 
Under Assumption 2 (which allows us to choose a contin- 
uously differentiable without loss of generality), Lemma 
14, and the facts that, |w(a;,t)| < J Q \w x (x, t)\dx, for 
all x G [0,1] (which follows from (31)), and that / : 
C 2 (R n x M; M n ), /(0, 0) = 0, which allows us to conclude 

\f{X,u)\<a 6 (\X\ + \u\), (B.32) 

for some function as € /Coo H C 1 , we get from (A. 2) that 



r(x,t)\ < -jrlukfot)! +a r {A{t),R) 



D 



(B.33) 
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A(t) = \X(t)\ + [ a*(\w(x,t)\)dx 
Jo 

- / \w x (x,t)\dx, (B.34) 
Jo 



+ 



for some class ICC oo function a r , continuously diffcrcn- 
tiable in its first argument. Analogously, differentiating 
(A. 2) with respect to x and using (28) together with the 
fact that Dr(x, t) = u x (x, t) it is shown that 

\r x (x,t)\ < jr\w xx (x,t)\ +a ltrx (A(t),R) 

+ (fJ-*+a 2 , rx {A(t),R))\r(x,t)\, (B.35) 

for some positive constant \x* and some functions ai ;rx , 
0-i,r x € K,Coo which are continuously differentiablc 
with respect to their first argument. With the Cauchy- 
Schwarz inequality we get (B.25), (B.31). 



C Proof of Lemma 2 

The proof of this lemma is based on the following fact. 

Fact 1 There exists a class AXoo function £i such that 
if the perturbation 8 and the solutions of the system sat- 
isfying (17) for < c < 1 and (50), the following holds 
for all 4>{t) <9 <t 



\P*(0)\<Ci [\X(t)\+ sup \U(s)\,R). (C.l) 

y 4>{t)<s<t ) 

Proof The proof can be found in [5] (Lemma 4). 

It holds U{6) = U(t) - f* U(s)ds, for all 4>{t) < 9 < t, 
and hence, using (20), (2) and the fact that y e U{s)ds = 
flw-t u x (x,t)dx we get sup 0(t)<e<t \U{6)\ <a (\p(l,t)\) 

+ J Q (\u x (x, t)\ + \u x (x, t)\) dx. Lemma 12 (App. B) and 
the Cauchy-Schwarz give 

sup \U(9)\<a 4 (\X(t)\+ [ a*{\u(x,t)\)dx 
<j>(t)<e<t V Jo 

ii x (x,t) 2 dx+ u x {x, t) 2 dx I (C.2) 
o Jo J 

for some class /Coo function a 4 . Using (22), (B.32), con- 
ditions (17) for < c < 1 and (50) are satisfied, if the 
following relation is satisfied for all cf){t) < 9 < t 



R 1 >c 1 + (i{\P'{6)\) + (ci + M\P'W)\)) 
xa 5 l\X(t)\ + sup |*7(0)|), 



(C.3) 



where i?i = min |c, Z)j. With Fact 1, Lemma 5 and 
(C.2) the lemma is proved with c* = c* 2 and a satisfying 



R> ( Ci + a(6 («4 {3( 2 {c*2,R)),R))) 
x(l + a 5 ( ai (3( 2 (c* 2 ,R))))- 

D Proof of Lemma 3 



(C.4) 



Let c* be the minimum of c\ and c\ defined in (B.20) 
and (C.4) respectively. Taking the derivative of V with 
S* = V~S and using integration by parts together with 
(29)-(30), (34)-(35), (36)-(37), (39)-(42), (B.18) we get 



+ 



M 4 L*(R) 



u(0,t)|-$ii7r(0,t)|u(0,t)| 



-giQii I e 9lX ir(x,t)\u{x,t)\dx - guir x (x,t) 
Jo 

: / e 9lX \u(x,t)\dx + gn sup 1 - L>TT(x,t) 
Jo xe[o,i] 



x / e 9lX \r{x,t)\dx + g 6 e 92 



tt(M) 

xr(l,t) 2 - 56 ^(0, t)ti x (0,t) 2 - g 2 g 6 



D 



x / e 92X ir(x,t)u x (x,t) dx + g 6 TT x (x, t) 
Jo 



f 

Jo 



x / e 92X u x (x,t) 2 dx + 2g 6 sup 1-Dir{x,t) 



f 

Jo 



ce[o,i] 



x / e 92X \r x {x,t)\\u x (x,t)\dx - g 7 w x (0,ty 



-9497 / e 9iX w x {x,t) dx-g 9 g w 
Jo 



f 

Jo 



x / e 9loX \w x (x,t)\dx-g 9 \w x {0,t)\+2g 7 



x / e 9 * x \w x (x,t)\\r 2 (x,t)\dx f(t) 



+99 



Jo 



r 2 (x,t)\dx\f(t)\+g 7 e 9 * ri (l,t)' 



f(t) +g 9 e^\n(l,t)\\f(t)\+g s e^ 

xw xx (l,t) 2 - g 8 w xx (0,t) 2 - g 8 g 5 [ e 95X 

Jo 

xw xx (x,t) 2 dx + 2g 8 / e 95X \r 3 (x,t)\ 
Jo 

x\w xx (x,t)\dx f(t) -g 12 a* (\w(0,t)\) 

-91293 a* (\w(x,t)\) dx + g 12 
Jo 
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f e 9sX \a*' 
Jo 



(\w(x,t)\)\\n(x,t)\dx f(t) , (D.l) 



for an increasing function L* £ C (R+;R+). Using (38) 
and Lemma 2 we get for all x £ [0,1], -. 4-* r < 

(l-\-R) ( D-\-R J 

^ { i-R)\b- R y M*,*)! ^ {1 -H)\b-Ry where 

c = (ci + A (Ci (o4 (3C 2 (min {cj, c^} , R)) , i?)) ) 

x(l + a 5 (a 4 (3C2(niin{ C t,^}, J R)))). (D.2) 

Moreover, since n(x, t) is linear in x, it takes its maxi- 
mum value either at x = or at x = 1 , and hence, 

l-Dir{x,t) < max j |l-£>7r(0,i) , 1 - £>tt(1, i) 1 1 
<M 2 Z(t), (D.3) 



where M 2 = 7r(0,t) + 7r(l,t) < ^Jg^v , and 

Z{t) = max{\5(a(t),P*(t))\, 8 {a{t), P* (t)) + D 

x(\5 t (a(t),P*(t))\ + \VH<r(t),P*(t))f(P*(t),U(t))\)} 
Therefore, using (C.4) we have that 



sup 

x£[0,l] 



1-Dir(x,t) <2cB*(R) 



B*(R) 



D + 2 



D-R l-R 



(D.4) 
(D.5) 



Since / : C 2 (R n x R; MP), using relations (32), (43) and 
(44), with Lemma 3.1 from [26] and (51) we have that 



f(t) + f p (t) + U(t) < Kl (R)\u(0,t)\, 



(D.6) 



for an increasing function n\ £ C (R+; R+). Hence, from 
(40), (B.33) we conclude after using (D.6), (B.34) that 



r(l, tf < k 2 (R)\u(0, t)\ + 2a 2 r (A(t),R) , 



(D.7) 



for an increasing function k 2 £ C (R+;R+) (where we 

2 



also used the fact that 



which follows form 



< c(R) 

(32) and (51)). We are concerned next with w xx (l,t) 2 . 
With Young's inequality and (B.25), from (42) we get 
that there exists an increasing function n 3 £ C (R+; R+) 



such that (where we absorb the powers of 



and 



higher than one in k 3 based on (32), (43), (44), (51)) 



fu 



w xx (l,t) 2 <n 3 (R)(f(t) + fp(t) + U(t) 
+k 3 (R)u x (0, tf + k 3 (R)w x (0, t) 

+ ( sup 1-Dir{x,t) I k 3 (R) 
\xe[o,i] / 

xa 2 (A(t),R) , 



(D. 



where we used (B.34) and the fact that r(0,i) 2 < 
■j^w x (0,t) 2 + 2a 2 {A(t),R) which follows from (B.25). 
From (B.34), (B.35) and (51) we get 

[ r x (x, t) 2 dx < j w xx (x 7 t) 2 dx + k 4 (_R) 
Jo D 2 Jo 

w x (x, t) 2 dx + a 2 r (A(t), R) 



+3a 2 (A(t),R). 



(D.9) 



With relation (51) and Lemmas 4, 5, 16 (App. B), from 
(D.l) one can conclude that the terms that multiply |/| 
are bounded by {g 7 e 91 + gge 910 + g$e 95 + gi2& 93 ) Ks(i?). 



Choosing. 9l =.g 2 = (l + i?) (b + R ) (l + 

and combining (D.7), (D.8), (D.9), (D.6) with Young's 
inequality, we get from (D.l), (B.33), (B.34) and (B.35) 



v(t)< 



' /: |Jr(t)| + ^^(|»a(o,t)| + |fi(o,t)|) 



2VMT 



-giij \u{x,t)\dx - gi 2 a* (|w(0,t)|) 
Jo 

-«/n(l + R)' 1 (b + Rj V(0, t)\ - u x (0, t) 2 

x (g^l + RY 1 (b + r)' 1 - g & e^n 3 {R)^ 

-(g 7 -g8e^ K3 {R))w x {0 7 t) 2 -g 6 

x (l-2cB*(R)) [ e 92X u x (x,t) 2 dx 
Jo 

-5i2ff3 / a* (\w(x, t)\) dx - g$w xx (0, t) 2 
Jo 

- (g 79i - 2ce 92 g 6 (bn A {R) + lJb~ 1 B*(R)^j 
x J e 9iX w x (x,t) 2 dx - (. 98 . 9 5 - g^D- 2 

x6cB*(R)) ( w xx (x,t) 2 dx 
Jo 

+ (^g 6 e 92 (^ + l) 2 |l- J D^(l,t)| 2 K2 (i?) 

+ {g 7 e 9i + g 9 e 910 + g 8 e 9s + ffl2 e 93 ) K 5 (R)) 
x|«(0,t)|- ^g 9 g w -giie 9 '2cB*(R)b-^ 

x f e 9wX \w x {x,t)\dx + 2g 6 e 92 cB*(R) 
Jo 



x (3a 2 !rx (A(t),R) + n A {R)a 2 r (A(t),R)) 

-g 9 \w x (0, t)\ +gu sup l-bn(x,t) 
xe[o,i] 



c9l 



xa r (A(t),R) + sup 1-Dn(x,t) (8g 6 
\xe[o,i] j 

xe 92 (b + l) 2 +g s e 95 K 3 (R))a 2 r {A(t),R) . (D.10) 
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From the proof of Lemma 16 (App. B) we have 
that a r and ai, rx are continuously differentiable 
and hence locally Lipschitz with respect to their 
first argument. Using (51) we can write a 2 (s,R) < 
a2e(R)a r (s,R) < L\(R)a2e(R)s and a\ Tx {s,R) < 
a2e(R)a r (s, R) < L2{R)a 2 &{R)s for every bounded s, 
some increasing functions Li <E C (K + ; K+), i = 1,2, 
and some class /Coo function a26- Choosing g 6 > 

(1 + R) (p + R)g 6 e!»K 3 (R), g 7 > g 8 e 95 n 3 (R), 9ll 



> 



(1 + R) i^D + Rj \4g 6 eP + 1 
+ (.97e 94 + g 9 e 910 + 9se g5 + g 12 e 



4B* 2 (R)n 2 (R) 
'<■■ i K 5 (R) 



MjL'(R) 
2VMT 



ff4 > £ (e 92 ff 6 («4(i2) + i) 2B*(i?) + 1), fl3 = ff5 = 
5 = 1 and since from the proof of 



9s = 59 



5io 



Lemma 14 a*(s) > s choosing gi 2 > ^t^i^ we get 



V(t)< 



^=-cB)\X(t)\- gil J o \u(x,t)\ 



u x {x, t) 2 dx 



xdx- g 6 (l-cB 2 ) f e 92X 
Jo 

- f w x {x,tfdx-{l-c{B 3 + B)) 
Jo 

x / \w x (x,t)\dx - (.9i2 - cB) 
Jo 

x a* (\w{x,t)\)dx - {I- cBx) 
Jo 

x / w xx (x,t) 2 dx, 
Jo 



(D.ll) 



where we used (B.34) and 

B (R) = 2g 6 e a2 B*{R)a 2(i {R) {3L 2 {R) + n^L^Rj) 
+32g 6 e g2 (f> + l) ii(i?)a 26 (i?)4B* 2 (R) 
2L 1 {R)B*{R) (g lie ai +2a 2(i (R)eB*(R))(D.12) 

B 1 {R)=6g 6 e 92 I)- 2 B*(R) (D.13) 



B 2 {R) = 2B*(R) 

B 3 {R) = 2eg 11 e 9l D- 1 B*(R). 



(D.14) 
(D.15) 



Restricting c* = min {c\, c^} and c\ such that c in (D.2) 
satisfies c<min {R, ci}, with cimax{_B2, B\,B 3 + B} < 



±min{ 



912, 



, l|, we arrive at V(t) < — AV(i), with 
{^cFT' 2 5n, 36, 1,312}- With the compari- 



2 mm 1 2VMI 



son principle (Lemma 3.4 in [26]) we get (55). 



E Proof of Lemma 4 



Using (B.14), Lemma 13 (App. B), the fact that for all 
x G [0, 1], |u(x,t)| < J Q \u x (x, t)\ dx (which follows from 



(35)), the Cauchy-Schwarz inequality and some routine 
class K- calculations the proof is immediate. 

F Proof of Lemma 7 

Taking the derivative of V L along (79)-(81), (83)-(84) 
and (89)-(90) and using integration by parts, we get that 

VL (t) < - \X(t)\ 2 A min (Q) + 2X(t) T PBw(0, t) 

+2X(t) T PBu(0,t)-b 1 b [ e bx n{x,t)u(x,t) 2 

Jo 

xdx- bin(0,t)u(0,t) 2 - h / e bx ir x (x, t) 

Jo 

xu(x, tfdx - 26i J e bx u{x, t) (l - Dir(x, tfj 

xr(x, t)dx + b 2 (^-w(0, t) 2 - J w(x, t) 2 dx 
+2D 2 \K\ 2 e 2 ^ i) x |B| 2 «(0,t) 2 - ^(0,t) 2 ) 

-b 2 [ w x (x,tfdx -2b 2 D [ (1 + x) 
Jo Jo 

xw(x, t)Ke Af)x Bu(0, t)dx - 2b 2 

x [ (l+x)w x (x,t)D 2 Ke Ai)x ABu(0,t)dx.(F.l) 
Jo 

With similar calculations as in [30] and with Lemma 6 
by choosing b > (1 — 7r") max jl, pU- j we get that 



bw(x,t) +n x (x,t) > tt*;/3*, 



(F.2) 



where /3* = min {b - 1 + ttJ* , (6 + 1)^* - 1} > 0. Since 
w(x,t) is linear in x, it takes its maximum value either 
at x = or at x = 1, and hence, 

(l - £>tt(M)) < max j 1 1 - -Dtt(0, t) , 1--ZMM)} 
< M 2 , L7 (t), (F.3) 



where 



M 2 ,l = 7r(0,t)+7r(l,t) < 7 



1 + sup f > to a(t) 
inf t > to (cr(t) - t) ' 



(F.4) 



and 7 is defined in (100). We derive next a bound for 
r(x, t) in terms of X, w and w x . Using (88) together with 
Young's and Cauchy-Schwarz 's inequalities we get that 

||rW|| 2 < M llL (\X{t)\ 2 + ||^(t)|| 2 + ||^(t)|| 2 ) (F.5) 



\W)f 



r(x, t) 2 dx, 



(F.6) 



Mi iL = 4I)- 1 + 41? Xe( A+SK ) £) (A + M) 



+ 4D 



x|^ J B| 2 +4Z) 



Ke( A+BK ^(A+BK)DB .(F.7) 
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Using (F.2), (F.3), (F.5) and Young's inequality we get 

VL(t) < - Ami " (Q) \X(t)\ 2 - [ e bx u(x,t) 2 dx 

4 Jo 



b f b f 

— \ w(x,t) 2 dx- — / w x (x,t) 2 dx 
2 Jo ^ Jo 

-b 2 )w(0,t) 2 +(^ 



8|PP| 



gjPgj 

-^min 

(Q) 

+86 2 L> 2 |X|el A l £, |B| 2 (2 + \A\ 2 ) - 6 l7 r**) 
x«(0, tf + 6iM 2 , L (l + Mi, L )7(t)H(t), (F.8) 



where 



(F.9) 



S(t) = |X(i)| 2 + / e bx u(x,t) 2 dx 
Jo 

+ / w(x,t) 2 G?a;+ / u> x (a;,f)' 
Jo Jo 

Choosing h = + 8b 2 D 2 \K\e\ A \ £) \B\ 2 (2 + \A\ 2 ), 

b 2 = x^f(Q) i anc ^ using the fact that 



M 4 , L H(t) < v L (t) < M 3 , L s(i), 

where 

max 

(P) + h + 2L>6 2 
M 4 , L = min {A min (P), 6i, £)6 2 } , 

we get relation (99) with 

= min{A min (Q),46i7r5*/3*,26 2 } 
ri 4M 3 , L 
fciM 2 x(l+Mi [L ) 

M 4 . L 



T 2 



(F.10) 

(F.ll) 
(F.12) 

(F.13) 
(F.14) 



T(0,t) 



. Using (87), we 



G Proof of Lemma 8 

Consider first that j(t) = 

§ et = 1h d+jSS'' ' and hence ' 7W = |J {am - 

Therefore, if 5 satisfies (65) with <5ir 2 < r\, the lemma is 



proved. Assume next ^(t) 



D 



Thus, 7(f) 



S{a(t)){l-S'(a(t)))-DS'(a(t)) . With (96), (65) we get 

j(t) < Si (l + Si + D^j . Hence, the lemma is proved if 

S\r 2 < tt, where d = 1 + Si + D. Note that since one 
can choose 7Tq* j3* sufficiently large (by choosing a large 
6) such that n in (F.13) is independent of (5i (or very 
large), and since from (F.14) together with (F.7), (F.4), 
(F.12) r 2 is bounded, one can always find a sufficiently 
small Si such that relation Sir 2 < is satisfied. 



H Proof of Lemma 9 

Using relations (77), (78) together with Young's and 
Cauchy-Schwarz's inequalities we get that 



H;)ll 2 <^i(l*(t)l 2 + IH*)!I 2 ) 
II^WII 2 < ^ (\x(t)\ 2 + \\w(t)\\ 2 + \\w x (t)\\ 2 ) 

\\w(t)\\ 2 <N 3 (\X(t)\ 2 + \\u(t)\\ 2 ) 

||^W|| 2 <^4(|A(t)| 2 + ||,i(t)|| 2 + ||^(t)|| 2 ), 

where || • || is defined in (F.6) and 



JVi = 

N 2 : 

N 3 -. 

N 4 : 



3 + 3 

:4 + 4 

+4L> 2 
3 + 3 

:4 + 4 

+4L> 2 



Ke (A+BK)b 



3L> 2 



Ke {A+BK ^ D {A + BK)D 



Ke (A+BK)D B 

2 



(H.1) 
(H.2) 
(H.3) 
(H.4) 



(H.5) 



Ke (A+BK)D 
2 



AD 



Ke 



Ke AD AD 



Ke AD B 



3D Z 



+ AD 2 \KB\ 2 
B\ 2 (l + D\A + BK\yR.6) 

Ke Ab B 2 



+ AD 2 \KB\ 



+ 4D 2 



Ke AD ADB 



(H.7) 



(H.8) 



Using (F.10), (H.1)-(H.4) the lemma is proved with 
M ^ = Tkz ( 2eb + N 3 + N 4 + 1), M 6 , L = ™±^. 
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